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Abstract 



We study the effect of confinement on the phonon properties of ultra-narrow silicon 
nanowires of side sizes of 1-lOnm Wg use the modified valence force field method to 
compute the phononic dispersion, and extract the density of states, the transmission 
function, the sound velocity, the ballistic thermal conductance and boundary scattering- 
limited diffusive thermal conductivity. We find that the phononic dispersion and the 
ballistic thermal conductance are functions of the geometrical features of the structures, 
i.e. the transport orientation and confinement dimension. The phonon group velocity and 
thermal conductance can vary by a factor of two depending on the geometrical features of 
the channel. The <110> nanowire has the highest group velocity and thermal 
conductance, whereas the <111> the lowest. The <111> channel is thus the most suitable 
orientation for thermoelectric devices based on Si nanowires since it also has a large 
power factor. Our findings could be useful in the thermal transport design of silicon- 
based devices for thermoelectric and thermal management applications. 



Index terms: Confined phonons, silicon nanowires, lattice thermal conductance, 
modified valence force field method, Landauer formula. 



I. Introduction 

Low dimensional silicon nanowires (NWs), and silicon-based ultra-thin layers 
have attracted significant attention as efficient thermoelectric (TE) materials after it was 
demonstrated that the thermal conductivity in such materials can be drastically reduced 
compared to bulk Si. Values as low as 7^ = 2W/mK (compared to 7f = 140W/mK for 
bulk Si) were achieved [1, 2]. This resulted in zr values close to zr~0.5, a large 
improvement compared to zr of bulk Si, ZT^^,;^ ~ 0.01 [1, 2]. The large reduction in the 
thermal conductivity was attributed to enhanced scattering of phonons on the surfaces of 
the nano channels. 

Numerous studies can be found in the literature regarding the thermal 
conductivity of Si NWs [3, 4, 5, 6, 7]. The effects of different scattering mechanisms, i.e. 
surface roughness scattering, mass doping, phonon-phonon scattering, and phonon- 
electron scattering have been investigated by several authors [8, 9, 10, 11]. In these 
works, it is demonstrated that the thermal conductivity in ultra-narrow Si NWs drastically 
degrades once the diameter of the NW is reduced below 50nm, or when scattering centers 
are incorporated. For even smaller NW diameters, i.e. below lOnm, the effect of 
confinement could further change the phonon spectrum significantly, and provide an 
additional mechanism in the reduction of the thermal conductivity [12]. This could 
provide additional benefits to the thermoelectric figure of merit ZT . In this work, we 
employ the modified valence force field (MVFF) method [13] to address the effect of 
structural confinement on the phonon dispersion, group velocity and ballistic thermal 
conductance of ultra-narrow Si NWs of diameters below lOnm. We consider different 
transport orientations, and different cross sectional sizes. We find that although the 
density of phonon modes does not show any orientation dependence, the phonon group 
velocities and sound velocities are strongly anisotropic. This results in a higher thermal 
transmission in the case of the <110> NW, which results in a higher ballistic thermal 
conductance compared to differently oriented NWs as well. On the other hand, the <111> 
oriented NW has the lowest thermal conductance, almost 2x lower than that of the <110> 
NW, and could be the optimal choice for NW TE devices. Our results, therefore, indicate 



that the proper choice of the NW geometry can provide a further reduction in the thermal 
conductivity of NW TE devices and consequently higher Z7 values. 



11. Approach 

For the calculation of the phononic bandstructure we employ the modified 
valence force filed method [13]. In this method the interatomic potential is modeled by 
the following bond deformations: bond stretching, bond bending, cross bond stretching, 
cross bond bending stretching, and coplanar bond bending interactions [13]. The model 
accurately captures the bulk Si phonon spectrum as well as the effects of confinement 
[14]. 

In the MVFF method, the total potential energy of the system is defined as [14]: 






(1) 



where JV^ , nn- , and COP- are the number of atoms in the system, the number of the 
nearest neighbors of a specific atom 7', and the coplanar atom groups for atom '/', 
respectively. L/^, , U^^, L/^..^,, (7^^.^^, and U^^_^^ are the bond stretching, bond bending, cross 
bond stretching, cross bond bending stretching, and coplanar bond bending interactions, 
respectively [14]. 



In this formalism we assume that the total potential energy is zero when all the 
atoms are located in their equilibrium position. Under the harmonic approximation, the 
motion of atoms can be described by a dynamic matrix as [15]: 



d = [dL]^ 




(2) 



where dynamic matrix component between atoms 7' and 'j' is given as shown in Ref. [14] 
by: 
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is the second derivative of the potential energy after atoms '/' and 'j' are slightly displaced 
along the m-axis and the n-axis (9rJ, and dr^ ), respectively. To calculate L/^,^^;;^, we added 
only those terms of Eq. 1 that include atoms 7' and 'j'. All other terms cancel each other 
out when we numerically calculate the second derivative of potential energy. 



After setting up the dynamic matrix, the following eigenvalue problem is solved 
to calculate the phononic dispersion: 

D + Y^D^ exp(/q.A^)- o/ (q)/ = 

(5) 
where D, is the dynamic matrix representing the interaction between the unit cell and its 

neighboring unit cells separated by ^R, [15]. Using the phononic dispersion, the density 
of states (DOS) and the ballistic transmission (number of modes at given energy) are 
calculated as shown in Ref. [16] by: 

DOS{o))=Y^DOSioj)=Y^d{o)-o,,{q)) 

(6) 
and 
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Once the transmission is obtained, the ballistic lattice thermal conductance is 
calculated within the framework of the Landauer theory as [17]: 
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It is well known that the thermal conductivity in Si nanowires of very narrow 
cross sections is dominated by boundary scattering, which is the reason of their improved 
thermoelectric performance as well [2, 11, 18, 19]. The diffusive thermal conductivity 
can be calculated using the phonon lifetime approximation in the phononic Boltzmann 
transport equations as [18]: 
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where v^.(q) is the group velocity of phono ns of wavevector q in subband / , given by 
v^.{q) = dct>.{q)/dq , and z",((7) is the scattering time. Here we compute the boundary- 
scattering-limited thermal conductivity by evaluating the boundary scattering time as: 

1 ri-pVg,,(q) 
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where d is the effective diameter (defined as d = yfiw where w is the side length of the 
NW cross section) of the NWs and P is the specularity parameter that takes values form 
to 1, determined by the details of the surface and surface roughness [18, 20]. For a 
smooth surface P = 1 , the scattering time is infinite (meaning zero scattering rates on the 
surface), and boundary scattering is fully specular. On the other hand, for a very rough 
surface P = , which results in fully diffusive boundary scattering. Recently, studies 
showed that experimental results for boundary scattering in narrow NWs are only 
explained once an almost fully diffusive boundary is assumed [21, 22]. 



III. Results and Discussion 

In this work we examine the effects of transport orientation and confinement 
length scale on the thermal properties of Si-NWs with square cross sections {W = H) and 
side sizes up to lOnm. We study the phononic dispersion, density of states, ballistic 
transmission, average group velocity, and lattice thermal conductivity of different 
structures. 



The phononic dispersions of Si NWs of 2nmx2nm cross section in different 
transport orientations, <100>, <110>, and <111>, are shown in Fig. 1. There are 
differences in the dispersions, especially in the low energy, low momentum region, which 
indicates that the thermal properties could be orientation dependent as well. Figure 2 
shows the ballistic transmission function of a different set of NWs, of side sizes 
W = H = 6nm in the same three orientations, calculated as indicated above in Eq. 7. The 
transmission function of the <110> NWs is the highest in most part of the energy 
spectrum, whereas the transmission of the <111> NWs is the lowest in almost the entire 
energy spectrum. As a result, the ballistic lattice thermal conductance of the NWs (see 
Fig. 3-a) calculated using the Landauer formula as specific in Eq. 8, shows that the 
<110> NWs has the highest thermal conductance compared to the NWs oriented along 
the <100> and <111> transport directions. The thermal conductance is the lowest for the 
<111> NWs of all the side sizes we have considered, all the way up to lOnm as shown in 
Fig. 3-a. The difference between the thermal conductance of the <110> NW which has 
the highest performance, and the <111> NW which has the lower performance, is a factor 
of ~2. Another observation is that the thermal conductance increases as the cross section 
of the NW increases. This is expected since the larger NWs contain more transport 
modes. The increase is close to linear. Once the conductances are normalized by the cross 
section area of the NWs, however, the resultant normalized conductances remain almost 
constant. This is indicated in the inset of Fig. 3-a. In this case, again, the <111> has 
clearly the lowest conductivity, almost 2 times lower than the <110> NW for all cross 
section sizes. In Fig. 3-b we show the boundary-scattering-limited lattice thermal 
conductivity for the nanowires of 2nm and lOnm cross section sides, versus the 
specularity parameter P. As P increases from (fully diffusive surfaces) to 0.5 (semi- 
diffusive surfaces), the conductivity naturally increases. The scattering is stronger for 
narrower NWs. For the NWs with cross section sizes of 2nm, thermal conductivities as 
low as ~2 W/mK are achieved for P=0, raising to ~5 W/mK for P-O.S. For NWs with 
cross section sizes of lOnm, the thermal conductivity is almost 5 times higher, at most 
~30 W/mK. These values are significantly lower than the thermal conductivity of bulk Si 
(l40W/mK), and once other scattering mechanisms are considered (i.e. phonon-phonon 
scattering and impurity scattering) they could lower even further. Interestingly, the order 



of the conductivities for NWs of different orientations is the same as in the ballistic case 
of Fig. 3-a, with the <110> NW having the highest and the <111> NW the lowest 
conductivity. The same anisotropy behavior is observed when considering only phonon- 
phonon scattering mechanisms [23], and different fully diffusive transport formalisms 
[18, 19], further demonstrating the importance of the phonon bandstructure in 
determining the thermal conductivity of Si NWs. For thermoelectric applications, 
therefore, where low thermal conductance is advantageous, the <111> NW would be the 
optimal choice as we show below. 

We have recently studied the role of transport orientations and diameter on the 
thermoelectric power factor of NWs using atomistic bandstructure simulations for 
electron transport [24]. With regards to anisotropy, in the case of n-type NWs only a 
small orientation dependence was observed. In the case of p-type Si NWs, however, we 
showed that the <111> NWs have significantly higher power factors compared to 
differently orientated NWs. The <111> Si NW channel is, therefore, the most 
advantageous for p-type channel thermoelectric applications since it provides 
simultaneously the highest power factor and lowest thermal conductivity compared to p- 
type NWs of different transport orientations. Using the thermoelectric power factor of 
Ref. [24] (Fig. 3c and Fig. 7c of Ref. [24]) and the calculated thermal conductivity, in 
Fig. 4a and 4b we show the ZT figure of merit for n-type and p-type NWs, respectively, 
versus the NW cross section size. For n-type NWs, the <111> direction is the most 
advantageous, providing ZT of ~0.5. On the other hand, a ZT value around 1 is calculated 
for a very narrow p-type <111> NWs as shown in Fig. 4b, almost 3 times higher than that 
of the <110> NW. This is partly due to its higher power factor originating from its large 
electronic mobility [25] and partly due to its lower lattice thermal conductivity. 

For the rest of the paper, we revisit the ballistic results, and attempt to elucidate 
the reasons behind the anisotropy of the thermal conductance and conductivity in terms 
of the NWs' phonon bandstructure. The ballistic transmission function is directly related 
to the product of the density of phonon states and phonon group velocity. Below we 
investigate each of these quantities and their geometry dependence in order to provide 



insight into the orientation dependence of the thermal conductance. For this, we consider 
NWs with W = H = 6nm in the three different orientations under investigation. Figure 5 
shows the phonon DOS, which is nearly the same for all orientations in the entire range 
of frequencies. This can be understood by considering the fact that the DOS is mostly 
related to the density of atoms in the cross section of the NW. The density of atoms is 
mostly a bulk property and does not change drastically in Si-based structures. Each atom 
has three degrees of freedom (3x3 elements in its Dynamic matrix component), and 
therefore provides three eigenmodes. Although for ultra-narrow structures i.e. of 
H-W-2nm there could be some differences in the DOS due to finite size effects a 
6nmx 6nm Structure is large enough such that the DOS is orientation independent. 

On the other hand, as shown in Fig. 6, the sound velocity, defined by the slope of 
the acoustic modes near the Brillouin zone center, is not an isotropic quantity. We note 
here that the 2 lowest modes in phonon dispersion of NWs are flexural modes. The third 
mode is the one used to calculate the transverse sound velocity, whereas the fourth mode 
gives the longitudinal sound velocity. Our results, as shown in Fig. 6, are in good 
agreement with previous studies on NW sound velocities [14, 26]. The <111> NW has 
the highest longitudinal sound velocity (green-solid), whereas the <100> NW the lowest 
sound velocity for all NW cross section sizes. The transverse velocities are lower, and in 
that case, the <100> NW has the highest velocity. The values calculated here are 
different, and in general lower than the bulk sound velocities (dots in the right side of Fig. 
6). The velocities in NWs are lower, also in agreement with other reports, since 
confinement flattens the phonon spectrum and reduces sound velocities [27]. 

At first, following a bulk way of analysis, it seems that according to the ordering 
in the sound velocity one would expect that the thermal conductance should be higher in 
the <111> NWs, followed by the <110> NWs and then by the <100> NWs. As we show 
in Fig. 3, however, the order is different, namely that the <110> NWs have the highest 
thermal conductivity, followed by the <100>, and then by the <111> NWs. The fact is 
that the phonon dispersion is modified in nanostructures so that the bulk definition of 
acoustic and optical modes is no longer the relevant quantity. The sound velocity 



determines the dispersion only in a small part of energy spectrum ( ~ 5meV of out of the 
total ~ 65meV). In addition, at room temperature, all the phonons in the entire energy 
spectrum contribute to the thermal conduction. Therefore, the sound velocity alone 
cannot describe the NW phonon spectrum, which is an indication that the bulk treatment 
is inefficient for understanding thermal transport in NWs. What needs to be taken into 
consideration in order to correctly interpret the results, therefore, is the group velocity of 
all phonon modes in all energies. For this, we define the effective/average group velocity 
as follows: 

h^^^h^J^Sic-cM) (11) 

•.1 

where the sum holds over all the modes. The velocity of a phonon is in general a function 
of its subband, its frequency, and its q-point. The quantity in Eq. 11 effectively 
suppresses the subband and q-point indices and provides a quantity that is only frequency 
(or energy) dependent. Similar ""effective" quantities are also used in actual thermal 
conductivity calculations with adequate accuracy in the results [8, 28], although in our 
actual calculations we utilize all the information of the phonon spectrum. 



In contrast to the sound velocity, the average weighted group velocity over 
different modes can provide a better and more correct insight with regards to phonon 
transport. Figure 7 shows that the average group velocity of phonons for NWs in the 
/110\ is the highest, or high enough compared to the other orientations, in the largest 

part of the energy spectrum. The explanation behind this lies in the shape of the 
dispersions shown in Fig. 1. In the case of the <111> NW, the modes look overall less 
dispersive, i.e. more "flat". This results in the lowest average group velocity. The modes 
of the <110> NW, on the other hand, are more dispersive, i.e. less "flat", which results in 
the highest group velocity. By considering that all NWs with the same cross sectional 
area have the same phonon DOS, the differences in the group velocity explain the 
orientation dependence observed in the transmission functions and in the thermal 
conductance shown in Fig. 2 and Fig. 3, respectively. Eventually the transmission 



function of the <111> NW is ~1.25x lower than that of the <100> NWs and almost 2x 
lower than that of the <110> NWs. 



IV. Conclusions 

In this work we study the thermal properties of ultra-narrow silicon nanowires 
using the atomistic modified valence force field method for the computation of the 
phonon bandstructure. We extract the thermal properties using the Landauer formalism 
under ballistic, and diffusive boundary-scattering-limited transport considerations. We 
address the effects of structural confinement on the phonon dispersion, the phononic 
density of states, the phononic transmission function, the sound velocity and the effective 
group velocity. Our results show that differently oriented NWs can have up to a factor of 
~2x difference in their effective group velocity, transmission function and finally ballistic 
and diffusive thermal conductance. The <110> oriented NW has the highest ballistic and 
diffusive thermal conductance, followed by the <100> and finally the <111> NW. By 
taking into account the thermoelectric power factor, we show that the <111> orientation 
is the optimal choice for NW TE channels, since it provides the lowest thermal 
conductance and highest power factors for both n-type and especially p-type channels. 
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Figure 1 caption: 

Phononic dispersions of Si NWs of square cross sections with W-H-2nm for a) <100>, 
b) <110>, and c) <111> transport orientations. 
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Figure 2: 
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Figure 2 caption: 

Transmission function versus energy for NWs in different transport orientations. <110> 
NWs have the highest and <111> NWs have the lowest transmission. 
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Figure 3: 
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Figure 3 caption: 

(a) Ballistic lattice thermal conductance versus NW cross section side size for NWs in 
different orientations. Inset: The thermal conductance normalized by the cross section 
area, (b) Diffusive boundary-scattering-limited lattice thermal conductivity for NWs with 
different orientations versus the specularity parameter. Nanowires with side sizes of 2nm 
and 10 nm are considered. 
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Figure 4: 
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Figure 4 caption: 

The thermoelectric figure of merit TT for (a) n-type and (b) p-type NWs with different 
orientations versus the nanowire side size. For this calculations, we used the 
thermoelectric power factors from Ref. [24]. 
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Figure 5: 
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Figure 5 caption: 

The phonon DOS versus frequency for NWs in different orientations. The DOS is nearly 
the same for the different transport orientations. 
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Figure 6: 
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Figure 6 caption: 

The sound velocity of NWs in different orientations versus the NW side size. 
Longitudinal and transverse velocities are shown in sohd and dashed lines, respectively. 
Dots: The respective bulk Si velocities. 
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Figure 7: 
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Figure 7 caption: 

Average group velocity of NWs in different transport orientations versus frequency. 
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